Statistical analysis of Ni nanowires breaking processes : a numerical simulation study 
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We have performed the statistical analysis of the breaking behavior of Ni nanowires. Using molec- 
ular dynamic simulations based on the embedded atom method approximation, we have simulated 
thousands of nanocontact breakages, determining the time evolution of both the nanowire atomic 
structure and the nanowire minimum cross section time evolution, S m (t). The accumulation of 
thousands of breaking events allows the construction of S m histograms. These histograms allow 
the analysis of the influence of the temperature, the crystalline stretching direction and the initial 
nanowire size on the breaking processes. We have also calculated the proportion of monomers, dimers 
and more complex structures at the latest stages of the breaking process. We have found important 
differences among results obtained for different nanowire orientations and sizes, illustrating that 
very small wires do not capture the behavior of actual thick nanonecks obtained in experiments. 
Three main cases have been observed: (a) For the [111] stretching direction and large nanowire 
sizes, the wire gradually evolves from more complex structures to monomers and, in several cases, 
to dimers prior its rupture. This case presents well ordered structures during the breaking process, 
(b) On the other hand, for large nanowires stretched along the [100] and [110] directions the system 
mainly breaks directly from complex structures, with a low probability of finding monomers and 
dimers. In this case a clear dependence on temperature is observed. At room temperature a huge 
histogram peak around S m — 5 appears, showing for first time in the literature the presence of long 
staggered pentagonal Ni wires with ...1-5-1-5-1-5-.... structure. Additionally, we have found for this 
case that nanowires have disordered regions during their breakage, (c) When the initial wire size 
is small strong size effects are observed independent on the temperature and stretching direction 
used. Finally, we have found that the local structure around monomers and dimmers do not depend 
on the stretching direction and that these configurations differ from those usually chosen in static 
studies of conductance. 
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I. INTRODUCTION 

Nanotcchnology involves the design, synthesis, fabri- 
cation, characterization, and application of structures 
by controlling their composition, shape and size at the 
nanometer scale. The control of these properties will 
allow to exploit a whole set of novel physical and chem- 
ical features in future devices with a wide range of ap- 
plications. In particular, electron transport through fu- 
ture nanofabricated devices is a crucial issue, and dif- 
ferent candidates have been proposed as potential nano- 
electronics elements with conducting or semiconducting 
properties. Among these candidates, metallic nanowires 
and nanocontacts are promising disipationless nanode- 
vices as well as exigent benchmarks to test nanoscale ba- 
sic phenomena 1 . 

The interest on metallic nanowires and nanocontacts 
rises from their very rich phenomenology. Electron trans- 
port through metallic nanowires present ballistic features 
at and below room temperature (RT) since, in general, 
the electron inelastic mean free path is larger than the 
characteristic nanowire dimensions. Furthermore, well 



defined electron transport modes or channels appear as- 
sociated to the transversal confinement of electrons for 
nanowires with diameters of the order of few Fermi wave- 
lengths Xp. In fact, a narrow metallic nanowire is a type 
of mesoscopic system that is characterized by the exis- 
tence of many electron levels (whose exact distribution 
depends on the precise atomic configuration) coupled to 
the electron energy spectra of the reservoirs. In the bal- 
listic limit, the conductance G (inverse of the nanowire 
resistance R, G = 1/R) is described in terms of the Lan- 
dauer formula 2 -: 



G=(G /2)J2J2 T - 



(1) 



where Go=2e 2 /h is the conductance quantum (e being 
the electron charge and h the Planck constant), a de- 
notes the two possible spin polarizations, N a is the spin- 
dependent number of accessible transport channels with 
energies below the Fermi energy, Ep, and T na is the spin- 
dependent transmission probability associated to the n- 
th mode. The calculation of the quantities T na is a diffi- 
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cult task since the solution to the scattering problem re- 
quires the full description of nanocontact band structure 
as well as the knowledge of those scattering mechanisms 
(geometry changes, defects, magnetic domain walls, etc) 
governing the electron transport. Consequently, the elec- 
tron transport is governed by the set of the non negligible 
transmission probabilities associated to the propagating 
modes. Notice that for degenerate energy levels Eq. (fTJ) 
is written as G — Gq~^2, x T n , where spin dependences 
disappear. 

Several experimental techniques have been used to 
form metallic nanocontacts and nanowires^. Scanning 
tunneling microscopy (STM)£4££ and mechanically con- 
trollable break junction (MCBJ) 7 -^ methods are stan- 
dard approaches to study, in a rather controlled way, 
the process of formation and rupture of nanocontacts un- 
der different experimental conditions. Metallic nanowires 
are also obtained using electron-beam irradiation inside 
ultra high vacuum (UHV) transmission electron micro- 
scopes (TEM) 9 .. This procedure allows the direct moni- 
toring of the metallic nanowires thinning process upon 
external perturbations (stretching, heating, etc). In 
some cases, this technique has revealed that nanowires 
with few atoms in the cross section present helical and 
weird structures^ -. Metallic nanowires are also formed 
through electrochemical methods by controlling the po- 
tential between pairs of electrodes immersed in elec- 
trolytic solution a 11 ' 12 . Finally, methods based on 'table- 
top' experiments (where contacts are broken by separat- 
ing two macroscopic metallic wires) 13 ' 14 provide a simple 
but uncontrolled way to generate metallic nanowires. 

The study of both ballistic and quantum features 
associated to the electronic transport through metal- 
lic nanowires requires their full electric characterization. 
This is usually done during the formation and/or the 
rupture of a metallic nanowire, allowing the access to 
many atomic configurations within the same experiment, 
and obtaining the nanowire conductance time evolution 
during its formation (or breaking) process, i.e. the so- 
called conductance trace, G(t). Each conductance trace 
presents its own features since it is very difficult to accu- 
rately control the nanowire geometry during its mechan- 
ical deformation. 

In order to obtain relevant information concerning 
the electronic transport through nanowires of a given 
metallic species, conductance histograms 6 - H(G) are con- 
structed by accumulating hundreds of conductance traces 
obtained under the same experimental conditions. Con- 
ductance histograms H{G) provide valuable statistical 
information on the transport and structural nanocon- 
tacts properties. In many cases, conductance histograms 
present well defined peaked structures close to integer 
multiples of Go , reflecting the existence of preferred 
conductance values. Such preferential values are usu- 
ally interpreted in terms of conductance quantization 15 
or favorable atomic arrangement a 16 ' 17 ' 18 . 

The interpretation of conductance histograms is a dif- 
ficult task since they merge mechanical and electrical 



information. Moreover, both properties are strongly 
correlated 19 . In fact, those atomic rearrangements taking 
place during the nanowire stretching process modify the 
number of accessible modes as well as the set of trans- 
mission coefficients, leading to conductance changes. The 
interpretation of conductance histograms is even more in- 
trincated for polyvalent metals since several channels per 
atom are involved in the electronic transpor t 20 ' 21 . For 
instance, aluminum conductance histograms obtained at 
both 4K 16 and RT 18 show well defined peaks at con- 
ductance values close to integer values of Go- How- 
ever, the first conductance histogram peak appearing at 
G ~ Go presents contributions from three propagating 
channels 20 . 

In spite of these interpretation problems, conductance 
histograms have become a standard tool to analyze the 
dependence of the nanowires electron transport on sev- 
eral experimental parameters as temperatur o 15 ' 22 , ap- 
plied voltage^ 3 -, adsorbed species^, and applied elec- 
trochemical potential 2 ^.. Furthermore, conductance his- 
tograms have been used to detect magic structures for 
different metallic nanowire o 15 ' 22 ' 26 ' 27 ' 28 . 

The simultaneous control of charge and spin 
currents using a suitable combination of geome- 
tries, applied magnetic fields and bias voltages 
opens an interesting landscape for future appli- 
cations. For this reason, electronic transport in 
atomic sized magnetic nanowires has been profusely 

^.j.^^5, 6 , 12, 29, 30. 31, 32, 33, 34, 35, 36, 37,38, 39, 40,41, 42,43,44, 45,46. 47,48, 49, 50, 51 , 5^ 

However, the analysis of conductance histograms H(G) 
becomes even more intricate in magnetic nanowires due 
to the presence of this new degree of freedom as well as 
the presence of new scattering sources as magnetization 
domain walls. 

There are substantive studies on nickel nanowires, in 
comparison to those referred to cobalt or iron, for in- 
stance. In particular, several experimental works have 
addressed the study of nickel conductance histograms 
under different environmental conditions. Although for 
most metals (notably Au) the histograms reported by dif- 
ferent groups are very similar, there exists a profusion of 
different H{G) curves for Ni. Therefore, it seems that Ni 
conductance histograms strongly dependen on the exper- 
imental conditions as well as on the nanowire fabrication 
methodology. 

The first study where a conductance histogram was 
reported 6 - also claimed the existence, at RT, of a Ni con- 
ductance histogram very similar to that of platinum, 
i.e. showing a peaked structure. These histograms 
were constructed using few conductance traces and ap- 
plying a selective criterion to choose them. However, 
further experiments (performed at RT, without applied 
magnetic field, and by adding thousands of unselected 
conductance traces) obtained featureless Ni conductance 
histogram a 30 ' 31 

Very different results, showing well defined peaks, were 
reported for Ni conductance histograms^ 2 - obtained in 
UHV, at room and higher temperature, and constructed 
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from many conductance traces. These histograms showed 
evidences of fractional conductance quantization as ex- 
pected for magnetic systems characterized by ballistic 
conductance and the presence on non-degenerated states 
(see Eq. [1]). More evidences of fractional conductance 
were obtained in several experiments carried out at RT 
where histograms were built up from few selected conduc- 
tance tracesi 33 ' 39 i 40 ' 49 Histogram peaks located at nonin- 
teger values of Go have also been observed during the 
formation of Ni nanowires by electrodepositioni^. Simi- 
lar results have been reported for other magnetic metallic 
specie o 34 ' 44 and for nanowires made of bulk non-magnetic 
metal s 14 ' 44 metallic nanowires. The later results could 
also point toward the presence of spontaneous magneti- 
zation in the nanowire 5 - 3 -. However, it has been recently 
demonstrated that the use of sets formed by few selected 
conductance traces allows the construction of histograms 
with fractional and non-fractional (i.e integer) conduc- 
tance peaks for Ni and Cu^. Therefore, conductance 
histograms formed with few conductance traces must be 
cautiously considered. 

A different set of experiments, carried out at 4K, 
in UHV conditions, and using thousand of conduc- 
tance traces, showed Ni conductance histograms with 
two well defined peaks around ~ I.6G0 and ~ 3.1Gq, 
respectively! 43 ' 46 ' 51 . The peaks were not modified by the 
presence of strong magnetic fields^ 6 -. The first peak po- 
sition is consistent with pioneering jump-to-tunnel ex- 
perimental results, obtained at T=4.2K— . More re- 
cent experiments 5 - 1 - show that the first peak is, in fact, 
formed by the superposition of two sub-peaks located at 
G - 1.2G and G - 1.5G . 

It is clear that there is a broad landscape with different 
Ni conductance histograms depending on the tempera- 
ture, the applied voltage, environmental conditions, etc. 
A similar situation happens for iron nanocontact o 31 ' 34 ' 37 . 
In summary, we can describe two opposite results for 
Ni when histograms are built without selecting conduc- 
tance traces: on one hand the low temperature conduc- 
tance histograms showing a well defined and reproducible 
peaked structure (not modified by the presence or ab- 
sence of applied magnetic fields), and, on the other hand, 
the room temperature conductance histograms, showing 
a bunch of different experimental results. The difference 
between low temperature (4K) and RT Ni conductance 
histograms is not well understood since, in both cases, the 
system is below its bulk Curie temperature (T C =903K) 
and a similar magnetic behavior is expected. 

The previously sketched experimental landscape con- 
trasts with more defined theoretical results. An infinite 
monoatomic Ni nanowire is characterized by a high con- 
ductance G = 3.5G(r^. This value decreases to G ~ Go 
in presence of magnetic walls due to the c£-states blocking 
effect^ 4 -. The same blocking effect is reported for more 
realistic Ni nanowires formed by three atoms attached to 
two (100) leads. In this case the conductance takes values 
comprised in the range 1.35Go — 1.6Go 55 >. It has been also 
reported that the conductance takes values in the range 



Go~2.5Go when different monomer and dimer configura- 
tions are proposed to describe the nanocontact region^ 9 .. 
More recent calculations on dimer-like Ni nanocontacts 
have demonstrated that G ~ I.8G0 for a full ferromag- 
netic situation whereas G ~ 1.4Go in presence of a do- 
main wall 5 ^. Therefore there are many theoretical results 
that could explain the first experimental peak appearing 
in Ni histograms at low temperature. However, the sec- 
ond peak of the conductance histograms remains unex- 
plained. 

Encouraged by the diverse experimental results, we 
have focused the present study on the role played by the 
mechanical behavior during the breaking process of Ni 
nanowires at low and room temperatures. Both tem- 
peratures are well below the Ni melting temperature 
(T m = 1728K) and, in principle, a similar mechanical 
behavior is expected. However, we need to exclude this 
thermal effect as the origin of those marked differences 
between low and room temperature H(G). The aim of 
the present work is to carry out a statistical study of 
the structural evolution of Ni nanowires under stretch- 
ing at low and room temperatures following a well tested 
methodology 17 ' 18 ' 27 ' 57 ' 58 . We will address the construc- 
tion of minimum cross-section histograms from hundreds 
of computational traces, taking into account different 
crystallographic orientations and nanowire sizes. We 
have not limited ourselves to the study of computational 
minimum cross section histograms but we have also an- 
alyzed the different types of structures appearing at the 
last stages of the breaking processes. An additional mo- 
tivation for carrying out the present study is also re- 
lated with the disagreement among different computa- 
tional histograms that have been recently reported for Ni 
nanocontact o 51 ' 58 ' 59 . These discrepancies must be under- 
stood in order to validate the computational construction 
of histograms as an useful methodology. 

The manuscript is organized as follows: in the follow- 
ing section we describe the computational approach used 
to construct computational histograms. In Section III, 
we present the dependence of such histograms on sev- 
eral parameters (temperature, nanowire stretching direc- 
tion, and nanowire size), paying attention to the differ- 
ent types of atomic configurations appearing at the last 
stages of the breaking process. Finally, the main conclu- 
sions will be described in Section IV. 



II. COMPUTATIONAL METHODS 

The study of those preferred configurations appearing 
in nanometric systems as well as their evolution under 
modification of external parameters is of capital inter- 
est. An accurate modeling of nanoscale systems will 
allow huge savings in expensive fabrication and char- 
acterization procedures. The main goal of this study 
is to understand the reation between the nanowire ge- 
ometry time evolution and the structure of the conduc- 
tance histograms. We follow a similar strategy to the 
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experimental one, simulating hundreds of independent 
breaking events, in order to determine the presence of 
preferred configurations leading to peaked structures in 
the conductance histogram. At this point, we should re- 
member that the quantity H(G)AG (where AG is the 
histogram bin width) is proportional to the time that a 
given conductance value G is experimentally observed, 
and, therefore, a conductance peak corresponds to a rel- 
atively stable conductance associated to one or several 
stable atomic configurations. 

The simulation of metallic nanowire breaking events 
has been carried out using standard Molecular Dynam- 
ics (MD) methods based on a semi-classical description 
of the atom-atom interaction. MD has been exten- 
sively used to study the structure and rupture of metallic 
nanowires. In fact, pioneering MD simulation s 60 ' 61 were 
used in the past to predict the abrupt changes in forces 
arising in a tip-surface indentation/stretching process. 

Ab-initio based methods have been used for describ- 
ing the dynamics of a single metallic nanowire under 
stretching 62 . In spite of their accuracy, the large com- 
putational resources required by ab-initio based methods 
make them unappropriated to determine computational 
histograms. A high computational effort is also required 
when other approaches are used as tight-binding molecu- 
lar dynamics (TB-MD) simulations 63 . It is also possible 
to determine computational conductance histograms tak- 
ing into consideration simple models^ In general, these 
methods do not account for a realistic treatment of the 
nanowire dynamics under stretching. Therefore, in or- 
der to perform the statistical analysis of many breaking 
events, semi-classical MD methods represent a midterm 
between accuracy and existing computational capabili- 
ties. 

In the past years, MD simulations dealing with metal- 
lic nanowires mainly focused on the description of sin- 
gle formation or breaking events, using different inter- 
atomic potentials but neglecting the study of statistical 
effecto 5 i 60 i 65 i 66 i 67 i 68 i 69 i 70 i 71 i 72 . More recently, several MD 
studies on breaking and formation processes of nanocon- 
tacts have statistically analyzed the appearance of pre- 
ferred atomic configurations in order to seek correla- 
tions with peaks noticed in experimental conductance 
histograms at different temperatures ! 17 ' 18 ' 27 ' 51 ' 57 ' 58 ' 59 ' 73 
Some of these MD simulation studie s 57 ' 59 ! 73 have been 
done using an hybrid scheme, where MD configurations 
are used as starting point for the calculation of the elec- 
tronic transport using accurate methods. These studies 
reveal that there exists a complex relation between the 
atomic configuration and the conductance histograms. 
For instance, for Al nanocontacts 5 - 7 - the first conductance 
peak around G ~ Go is caused by dimer-like configura- 
tions whereas the second peak close G ~ 2Go is due to 
monomer-like configurations. 

We have simulated the nanowire dynamics using a MD 
scheme where atomic interaction is represented by em- 
bedded atom method (EAM) potentials 7 ^. In EAM, the 



potential energy function for the system reads: 

E = \H + ( 2 ) 

where i and j run over the number of atoms. In the 
first term, <j)(r) corresponds to a pair potential depend- 
ing only on the distance between every pair of atoms 
i and j. The second term is the so-called embedding en- 
ergy, which depends on the mean electronic density pi 
at atom i's location. This density is approximated in 
EAM as the sum of the contributions due to the sur- 
rounding atoms, pi = J2j^i P( r ij)- Then, the embedding 
energy is calculated by evaluating and summing the em- 
bedding function F(p) at each atom's position. Depend- 
ing on the material and the specific physical properties 
to be studied, different pair potential 0(ry), embedding 
F(p) and density p{rij) functions can be defined. In the 
present study, we have used the EAM parameterization 
proposed by Mishin et ali&, which is constructed by fit- 
ting almost 30 experimental and ab-initio values corre- 
sponding to several bulk and surface properties. We have 
already used this EAM approach in previous works to 
describe static and dynamical properties of Al and Ni 
nanowiresj 17 ! 18 ' 27 ! 57 ' 58 ' 76 Although the simulation does 
not include explicitly the possible magnetization of Ni, 
its magnetic nature has been taken into account through 
the EAM potential since the embedding energy was fitted 
to some quantities obtained from spin-polarized ab-initio 
calculations^ 5 -. Furthermore, such magnetic effects are 
negligible on the nanowire dynamics due to the weakness 
of magnetic forces with respect to the crystalline field. 

Nanowire dynamics have been studied at constant tem- 
perature T using a standard velocities scaling algorithm 
at every MD step. Two temperatures (4 K and 300 K) 
have been considered in this work, in order to describe 
the Ni nanowire dynamics at low and room temperatures. 
The MD time step is St = 10~ 2 ps. Atomic trajectories 
and velocities were determined using conventional Ver- 
let velocity integration algorithms. The simulation of a 
single nanowire breaking event consists of three stages: 
i) the definition of a bulk-like unrelaxed nanowire; ii) a 
MD relaxation stage of the initial structure; and iii) a MD 
stretching process at constant velocity till the nanowire 
breaks. 

The first stage corresponds to the definition of the 
initial unrelaxed structure. We consider a bulk super- 
cell with parallelepiped shape that contains hundreds of 
atoms ordered according to a fee structure with bulk Ni 
lattice parameter a=3.52 A. The initial parallelepiped 
height (L) will correspond to the stretching direction and 
is larger than the base edges. We define the z axis as the 
stretching (pulling) direction. In the present study we 
have considered three different stretching directions cor- 
responding to the [111], [100] and [110] crystalline direc- 
tions. In addition, five different parallelepiped sizes have 
been considered for each stretching direction. Table U 
shows the geometric parameters characterizing the differ- 
ent initial bulk- like nanowires we have taken into account. 
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TABLE I: Geometric parameters of the initial non-stretched 
nanowires with parallelepiped shape. The total number of 
atoms in the nanowire (N a ) (including the frozen slabs) is 
calculated by multiplying the number of layers (Nl) which are 
perpendicular to the stretching direction [ijk] by the number 
of atoms per layer (N 3 ). 



For a given stretching orientation the three largest initial 
configurations have equal aspect ratio. The two small- 
est configurations have been chosen in such a way that 
all the configurations present similar initial cross section 
sizes independently on the stretching direction. The use 
of three different crystallographic stretching directions as 
well as five different initial cross-section sizes will allow 
to study the influence of these factors on the nanowire 
breaking mechanics at the two chosen temperatures. At 
the beginning of the simulation the velocity of each atom 
is assigned at random according to the Maxwellian dis- 
tribution that correspond to the simulation temperature. 

The second stage corresponds to the relaxation of the 
bulk-like initial structure. Firstly, we define two sup- 
porting bilayers at the top and bottom of the supercell. 
Notice that data shown in table |T] includes the supporting 
bilayers atoms. Atomic x and y coordinates within these 
bilayers will be kept frozen during the simulation. The 
nanowire will remain attached to these two bulk-like sup- 
porting bilayers during the relaxation stage. This stage 
lasts for 3000 MD steps in order to optimize the nanowire 
geometry. 

During the third stage, the z coordinate of those atoms 
forming the top (down) frozen bilayer are forced to in- 
crease (decrease) a quantity AL Z — 10~ 4 A after every 
MD step. This incremental process simulates the sepa- 
ration of the supporting bilayers in opposite directions 
(symmetric stretching) at a constant velocity of 2 m/s, 
giving rise to the subsequent nanowire fracture. Those 
atoms non located at the supporting slabs move follow- 
ing the EAM forces. Notice that the stretching velocity 
is much larger than that used in experiments. This draw- 
back is common for all the computational MD approaches 
(ab-initio or semiclassical). However, our computational 
description of the nanowire breaking event is compara- 
ble to the " slow" dynamics of actual experimental traces 
since the stretching velocity is smaller than the sound 
speed in nickel, allowing a nanowire evolution that fol- 
lows a set of metastable equilibrium configurations. 

During the stretching stage, the accurate knowledge of 
the atomic coordinates and velocities allows the deter- 
mination of many physical quantities characterizing the 



geometry and energetics of the evolving nanocontact. In 
particular, we have payed attention to the minimum cross 
section S m since it is closely related to the favorable con- 
figurations appearing during the nanowire evolution un- 
der stretching. Furthermore, S m is closely related to the 
maximum number of propagating channels and it deter- 
mines a first-order approximation of the conductance G 
using semiclassical approaches 7 - 7 -. 

The minimum cross section S m is calculated in units 
of atoms following standard procedures that have been 
successfully used in previous studie d 65 ' 66 . In our case, 
we define the radius ro to be equal to half the fee (111) 
interplanar distance (ro = dm/2). We assign a volume 
Vq = 47rrQ/3 to each atom. To calculate the section at 
a given z% position, we compute the total atomic volume 
Vtot,i inside a "detecting slab or cursor" with width Az. 
We have used Az = 4n- Notice that this cursor defini- 
tion has been also used in previous studies on nanowires 
stretched along different crystallographic direction o 59 ' 73 . 
In the present study we keep the same cursor size Az 
for all the nanowire crystalline directions. This allows a 
true comparison between histograms obtained for differ- 
ent orientations. The quantity Si — Vt t,i/Vo corresponds 
to the section (in number of atoms) at the position. 
The detecting cursor moves along the z-axis between the 
two frozen bilayers, using a step equal to 0.1 x dm. This 
scanning-like process determines the cross section Si as 
a function of z\ for a given computational time. Finally, 
from the set of collected Si values we determine the min- 
imum cross section S m at that time. In our study S m 
is calculated every 10 MD steps. We consider that the 
nanowire breaking process is completed when S m = 0. 

The evolution of S m with time shows a typical staircase 
trace as those examples shown in Fig. [T] These curves 
correspond to the largest simulated nanowires and they 
have been selected to illustrate different configurations 
appearing just before the nanowire breaking point. S m (t) 
traces show a stepped profile, with well marked jumps 
associated to atomic rearrangements of the nanowire. 
We have verified that these jumps are correlated with 
jumps in the force acting on the supporting slabs. In 
general, S m monotonically decreases between two sub- 
sequent jumps reflecting the existence of elastic stages. 
These elastic stages have been associated to the exper- 
imentally observed conductance plateaus. S m features 
larger fluctuations at room temperature than at 4K for 
the three studied nanowire orientations. 

The present computational method has been previ- 
ously tested to simulate nanowire rupture o 17 ' 18 ' 27 ' 57 ' 58 . 
However, this method poses some problems when used 
to detect the presence of atomic chains, dimers and 
monomers, since all these configuration are character- 
ized by S m values close to ~ 1. Moreover, the correla- 
tion between S m and conductance is unclear for low S m 
values^ 7 - ' 59 ' 73 . Therefore, in the present study we have 
focused our attention on the determination of S m as well 
as on the type of atomic configurations formed during 
the final steps of the breaking stage. We have found that 
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FIG. 1: Evolution of the minimum cross section S m for Ni 
nanowires stretched along the [111] (a-b), [100] (c-d) and [110] 
directions (e-f), at 4 K (a, c, e) and 300 K (b, d, f). This 
figure shows traces corresponding to the largest initial config- 
urations for each direction (1008, 1029 and 1015 atoms for the 
[111], [100] and [110] orentations, respectively). These S m (t) 
curves illustrate different breaking behaviours. Labels 'm', 
'd' and 'o' refer to monomer, dimer, and other structures; see 
text and Fig. [2] 



final stages can be classified according to three different 
categories: monomers, dimers and "others". These cat- 
egories are illustrated in Fig. where we have depicted 
the narrowest part of the nanocontact before its breakage 
for those curves (obtained at T=4K) shown in Fig. [TJ 

The monomer structure is characterized by a central 
atom standing between two "pyramids". This configu- 
ration usually shows a plateau around S m ~ 1. In the 
dimer structure, the apex atoms of two opposite pyra- 
midal configurations form a two atom chain. In some 
cases we have found that the dimer is not parallel to the 
stretching direction. This behavior is caused by the rela- 
tive shift between the two opposite pyramids during the 
nanocontact evolution. Depending on the relative angle 
formed between the dimer and the stretching direction 
we find values of S m ranging from ~ 0.5 to ~ 1.0. Final 
configurations, presenting more complex structures, that 
do not match these two categories have been labeled as 
"others". In this cases we generally observe an abrupt 
jump from S m > 2 (i.e., structures formed by two or 
more atoms) to S m = 0. As expected^, we have not 
found any atomic chain formed by three or more atoms 
for the three studied stretching directions. 

To detect the presence of monomers or dimers dur- 
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FIG. 2: Nanowire configurations just before breakage at 
T=4K and for three crystalline nanowire orientations [111], 
[100] and [110]. Structures were taken from those S m (t) 
traces shown in Fig. [T] Images illustrating the [111] orien- 
tation (a,b,c) correspond to times t = 0.77, 0.79, 0.74 ns for 
monomer, dimer and other configurations cases, respectively. 
For the [100] direction, snapshots (d,e,f) correspond to times 
t — 0.77, 0.83, 0.75 ns, respectively. For the [110] stretching 
direction, snapshots (g,h,i) were taken at t — 0.86, 0.87, 0.85 
ns, respectively. 



ing the ruptures we used the 'burning' algorithm intro- 
duced by Herrmann et al^S to study the internal struc- 
ture of percolating clusters on square lattices. This algo- 
rithm detects the presence of single connected bonds (also 
called 'red' bonds) in percolating clusters and networks. 
Monomers and dimers are composed by atoms that are 
equivalent to 'red' bonds if we consider the nanowire as a 
network of connected atoms. Classically, if we eliminated 
these particular atoms from the nanowire there will be 
no conduction between the external nanowire supporting 
slabs. We defined that two atoms are connected when the 
distance between their centers is lower that a quantity 
that we have called the conduction distance e?c- In par- 
ticular, in all results reported here we used dc = 3A We 
tested conduction distances closer to the nickel nearest 
neighbors distance but not relevant differences between 
results were observed. In our study we have applied the 
'burning' algorithm every 1000 MD steps in order to save 
computational time and, simultaneously, define a reason- 
able statistical ensemble. 
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FIG. 3: Minimum cross-section histograms H(S m ) (black 
filled curves) built with 300 independent Ni nanowire rup- 
tures under stretching along the [111] direction at T=4 K 
(a,c,e,g,i) and 300K (b,d,f,h,j). Different rows correspond to 
different initial parallepides sizes as indicated by the right 
side labels (atoms per layer x number of layers). For two 
nanowire sizes (1008 and 90 atoms) we also show the par- 
tial decomposition of the minimum cross-section histograms 
H(Sm)- Gray filled histograms were constructed taking into 
account those samples that showed monomers or dimers in 
the interval 0.25 < S m < 1.75. 



III. RESULTS AND DISCUSSION 

A. Minimum cross-section histograms 

Minimum cross section histograms H(S m ) have been 
built by accumulating S m (t) traces acquired during the 
simulation of hundreds of nanowire stretching and break- 
ing processes. In Figs. [3l [4] and [5] we show the histograms 
H(S m ), at 4 K and 300 K, for the [111], [100] and [110] 
stretching directions, respectively. Each histogram has 
been constructed with 300 independent breakages. The 
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FIG. 4: Minimum cross-section histograms H(S m ) (black 
filled curves) built with 300 independent Ni nanowire ruptures 
under stretching along the [100] direction at T=4 K (a,c,e,g,i) 
and 300K (b,d,f,h,j). Different rows correspond to different 
initial parallepides sizes as indicated by the right side labels 
(atoms per layer x number of layers). For two nanowire sizes 
(1029 and 90 atoms) we also show the partial decomposition 
of the minimum cross-section histograms H(S m )- Gray filled 
histograms were constructed as described in Fig. [3] 



different histograms are depicted in the range < S m < 
16. Note that this range includes the whole evolution 
of those S m (t) traces corresponding to the small size 
nanowires. These histograms clearly show the effect of 
temperature, initial nanowire size, and crystallographic 
stretching directions on the average structural evolution 
of breaking Ni nanowires. 

A first inspection of these figures reveals the ex- 
istence of well defined peaks associated to preferred 
nanowire configurations as it has been shown in previous 
work o 17 ' 18 i 27 ' 57 i 58 . A first general feature we can easily 
identify is the effect of temperature: peaked structures 
are rather sharp at T=4K whereas they present more 
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FIG. 5: Minimum cross-section histograms H(S m ) (black 
filled curves) built with 300 independent Ni nanowire ruptures 
under stretching along the [110] direction at T=4 K (a,c,e,g,i) 
and 300K (b,d,f,h,j). Different rows correspond to different 
initial parallepides sizes as indicated by the right side labels 
(atoms per layer x number of layers). For two nanowire sizes 
(1015 and 126 atoms) we also show the partial decomposition 
of the minimum cross-section histograms H(S m )- Gray filled 
histograms were constructed as described in Fig. [3] 



rounded shapes at T=300K. Small peaks at T=4K cor- 
respond to metastable configurations with slightly higher 
cohesive energies with respect to other metastable config- 
urations. The increase of temperature allows the explo- 
ration of more configurations during the stretching pro- 
cess, and, in this way, those metastable configurations 
with local minimum energy are easily accessible, leading 
to a better definition of their associated H{S m ) peaks. 

We have observed, in general, that the temperature 
increases preserves the general shape of Ni H{S m ) his- 
tograms leading to a better definition of the peaked struc- 
ture. This allows, for instance, an easier identification of 
those 'magic' configuration o 17 i 27 ' 58 . The thermal effects 
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TABLE II: Number of configurations (N con f) with minimum 
cross section S m in the interval 0.25 < Sm < 1.75, found in 
the corresponding set of 300 independent nanowire breaking 
events. N con j is shown for [111], [100] and [110] stretching 
directions, five nanowire sizes N a (see Table UJ) and two tem- 
peratures. 



are more remarkable for [100] and [110] stretching direc- 
tions where rounding effects of the peaked structure, and 
the disappearance of small peaks, become more evident. 
The more noticieble difference between low and room 
temperature histograms correspond to the [100] and [110] 
cases. For both orientations, histograms present an pro- 
truding S m = 5 peak at T=300K. We will discuss in the 
next subsection the structure of this peak which is not 
developed at low temperatures. 

A second streaking feature is that the overall his- 
togram shape is modified when the initial nanowire cross- 
section decreases. The first change is the obvious dis- 
appearance of high S m peaks as the initial nanowire 
cross-section takes smaller values. This happens because 
large S m values are not accessible for small initial sec- 
tion nanowires. The histogram H(S m ) also changes in 
the low S m region for decreasing initial cross-sections N s 
(see Table Q} . This change is more remarkable for the 
[111] histograms (see Fig. [3]) which are characterized by 
H(S m = 1) > H(S m = 2) > H(S m = 3) at T=4K for 
initial cross sections N s < 16. However this decreas- 
ing peaked structure changes at T=4K when N„ < 9. 
At higher temperatures (T=300K) the low cross-section 
region of the [111] histogram slightly changes when the 
initial nanowire thickness decreases, keeping the trend 
H(S m = 1) > H(S m = 2) > H(S m = 3). A similar situ- 
ation happens for [100] and [110] histograms: effects due 
to the initial nanowire sections are more evident at low 
temperatures. At room temperature H(S m ) histograms 
present similar peaked structures in the small S m region 

{S m < 2). 

Therefore, we have found that H(S m ) histograms are 
very dependent on the stretching (i.e the nanowire axis) 
direction. Whereas the [111] direction provides his- 
tograms with well defined decreasing peaked structure 
in the low S m region, the situation dramatically changes 
for [100] and [110] directions. These two orientations 
presents H(S m ) histograms with a clear 'depletion' in 
the region S m < 2. This is more evident when we ana- 
lyze the total number of breaking events (N con f) we have 
recorded in the interval 0.25 < S m < 1.75. This quantity 
corresponds to the area below the histograms curve in 



this S m range. The quantity N con f is shown on Table HT1 
for all the stretching directions, temperatures and initial 
nanowire geometry. 

It is worth see that the temperature determines two 
different behaviors of N con f as a function of the nanowire 
size. At low temperatures (T=4K), for large initial 
nanowire cross-sections, N con f is larger for the [111] 
nanowire directions when compared to the values found 
for the other two studied directions, [100] and [110]. How- 
ever, for the smallest initial nanowire cross-sections, the 
situation reverses and N con f is larger for [100] and [110] 
directions. At room temperatures (T=300K), the quan- 
tity N con f obtained for [111] nanowire directions takes 
larger values than those obtained for [100] and [110] 
nanowires in the whole range of initial nanowire sizes 
we have considered. 

The dependence of H(S m ) on the nanowire size and 
orientation is a very important issue since theoretical 
conductance histograms H(G) will also reflect this de- 
pendence. In our particular case (nickel), the use of 
the [100] stretching direction as well as narrow initial 
nanowires will determine a conductance histogram with 
small contributions from the low S m region (narrower 
nanocontacts), and in particular with low contribution of 
those peaks located within the region 0.25 < S m < 1.75 
(mainly formed by monomers and dimers as we will show 
later on this study). On the contrary, the [111] stretch- 
ing direction and a thick initial nanowire provide a dif- 
ferent conductance regime with large contributions com- 
ing from the region 0.25 < S m < 1.75. The important 
point we would like to remark is that these differences ex- 
plain the apparently conflicting histograms H(S m ) found 
in the literature. It is clear that calculated histograms 
corresponding to the [111] stretching direction, at low 
and room temperatures, and large initial sections^ (see 
Fig. E3a,b)) are rather different to that calculated for 
the [100] stretching direction, at T=4K, and small initial 
sections^ (see Fig. BJi)). 

Finally, we believe that, whatever it is the stretching 
direction, computational histograms obtained from large 
initial section nanowires must be taken into account in or- 
der to explain experimental results, where usually initial 
cross sections (i.e. conductance values) are very large. 
On the contrary, those histograms obtained from small 
initial cross section simulations should be carefully han- 
dled since they do not capture the behavior of actual 
thick nanonecks. 



B. Analyzing the peak at S m ~ 5 for [100] and [110] 
stretching directions: pentagonal Ni nanowires 

We have pointed out in the previous subsection that 
minimum cross-section histograms H(S m ) present a huge 
peak around S m = 5 for the [100] and [110] cases at 
T = 300K (see Figs. 0] and 03). If we restrict our study to 
thicker nanowires, it is clear that the largest well defined 
peak at S m = 5 corresponds to the [100] case. We have 
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FIG. 6: Distribution function P 4 . 8 _ 6 . 6 (At) at T = 300K of 
the time spent by the largest diameter nanowires during their 
breaking process within the minimum cross-section interval 
4.5 < S m < 5.5, for the three stretching directions: [111] 
(bottom), [100] (middle), and [110] (top). 



carefully analyzed the different nanowire configurations 
which contribute to the formation of this peak. First, we 
have calculated the distribution function P4. 5 „5. 5 (Ai) of 
the time At lasted by the nanowire between S m = 4.5 
and S m = 5.5 for the three stretching directions at T = 
300K. Results corresponding to P4.g-5.5(At) are shown 
in Fig. [6l It is clear that the distribution is wider for 
the [100] case, in correspondence with the information 
provided by H(S m ) histograms. Notice that we have 
found some S m (t) traces where the nanowires stays at 
the 4.5 < S m < 5.5 region for At ~ 200 ps. 

Our MD simulations allow to monitor the evolution 
of the full set of atomic coordinates during the break- 
ing process. In particular, we have depicted in Fig. [7] 
some examples of those nanowires with a relatively long 
staying-power time in the 4.5 < S m < 5.5 region, defin- 
ing the pronounced peak at S m = 5, for the [100], [110] 
and [111] cases. Note that the origin of the huge peak 
is related to the tendency of the system to form a long 
wire with S m ~ 5. At this point we should mention that 
these long structures are very common for [100] and [110] 
stretching directions, whereas rarely occur for the [111] 
case. In general, we have found that those nanowires 
forming long structures break without a progressive di- 
minishing of its atomic section, i.e. the system does not 
form monomers or dimers just before the rupture. On 
the contrary, when the nanowire forms a long structure, 
it usually breaks following the rupture pattern we have 
described as "others" (see Fig. [2j . 

A closer inspection to Fig. [7] allows to determine that 
these depicted wires present a well defined sequence of 

the atoms • 5 — 1 — 5 — 1 . This sequence does not 

correspond to any crystallographic fee or bec structure. 
This type of arrangement is not seen at 4K because a 
larger temperature is required to explore and overcome 
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FIG. 7: Snapshots of three configurations illustrating the ap- 
pearance of long wire structures with atomic section close to 
S m ~ 5 for the [110] (a), [100] (b), and [111] (c) stretching 
directions. These configurations correspond to the largest ini- 
tial nanowire cross-sections we have considered (see Table [TJ> . 
The rigth column shows cross-sectional perspective views as 
seen from the position indicated on the right column figures. 
The atomic configuration at the central part of the nanowire 
shows the staggered pentagonal structure -5-1-5-1- discussed 
on the text. 



those energy barriers leading to configurations able to 
develop these pentagonal chains. We have found that 
these pentagonal nanowires are formed by subsequent 
staggered parallel pentagonal rings (with a relative ro- 
tation of 7r/5) connected with single atoms. A similar, 
but shorter, structure was found for Cu(lll) breaking 
nanowires using MD simulations^. The stability of such 
pentagonal Cu nanowires was later confirmed by ab-initio 
calculations 80 , demonstrating that staggered pentagonal 
nanowires are favorable configurations. More recently, 
MD tight-binding calculations have determined the pres- 
ence of such pentagonal -5-1-5-1- nanowires during the 
breaking process^ of [110] Cu nanowires. Pentagonal 
motives also appear in infinite Al nanowires 8 -^ 

The most relevant point is that we have statistically 
demonstrated that the RT stretching of Ni nanowires 
along the [100] or [110] directions provide a method to 
generate staggered pentagonal nanowires with relative 
high probability. 



C. Analysis of the formation of monomer and 
dimer configurations 

The study of the minimum cross-section interval 
0.25 < S m < 1.75 is important because it defines the 
set of atomic configurations that present the type of 
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TABLE III: Total number of monomers (N mon ) and dimers 
(Ndim) analysed with the Herrmann 78 algorithm for the 
largest studied nanowires (with size N a ), and for [111], [100] 
and [110] stretching directions, at T=4K and T=300K. 



structures depicted in Fig. [2l i.e. monomers, dimers and 
more complex structures. The electron transport during 
the last nanowire breaking stages strongly depends on 
the type of nanowire atomic arrangements characterizing 
these final stages. Therefore it is critical to know the type 
of involved atomic configurations (monomers, dimers, 
other) appearing in the interval 0.25 < S m < 1.75 
(around the S m = 1 histogram peak) as well as their 
statistical weights. We have analyzed these statistical 
weights for three stretching directions ([111], [100] and 
[110]), two temperatures, and several nanowire sizes. 

The geometric characterization of the 0.25 < S m < 
1.75 region has been accomplished using the Herrmann 
'burning' algorithm every 1000 MD steps. Note that the 
number of configurations we analyze is smaller than the 
number of data used to construct the H(S m ) histogram. 
As an example, Table Hill shows the number of configura- 
tions where dimer or monomer have been detected for the 
three largest nanowires corresponding to [111], [100] and 
[110] stretching directions, respectively. A first inspec- 
tion of Table IIIII reveals that the number of situations 
with monomer and dimers is approximately five times 
larger for the [111] case than for the [100] and [110] cases. 
This means that, for the [111] case, the nanowire breaks 
gradually from more complex structures to monomers 
which, in some situations, form dimers. For the other 
two stretching directions there are many configurations 
that directly break from situations with minimum cross- 
section S m > 2. In consequence, the H(S m ) peak at 
S m = 1 is the largest one for the [111] case, whereas this 
peak is considerably smaller than those peaks formed at 
S m > 2 for the other two stretching directions. Interest- 
ingly, we have found for smaller initial system sizes that 
the number of configurations with monomers and dimers 
is similar for the three stretching directions. 

The fractions of monomers, dimers and complex struc- 
tures appearing in the range 0.25 < S m < 1.75 are de- 
picted in Fig.[H For the [111] stretching direction we note 
that the join proportion of monomer and dimers (~ 80%) 
is larger than the fraction of other structures (~ 20%). 
This fraction seems to be rather independent on the tem- 
perature as well as the initial nanowire size. This is a very 
interesting result since it implies that final configurations 
of a breaking Ni nanowire under [111] stretching mainly 
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FIG. 8: Fraction of monomer, dimers and more complex 
structures appearing during the Ni nanowire breaking pro- 
cess within the range 0.25 < S m < 1.75. 



correspond to monomer and dimer structures. 

For the other two stretching directions [100] and [110], 
the fraction of other (more complex) structures increases 
for both temperatures and also shows a larger depen- 
dence on the nanowire size. In some cases, the fraction 
of complex structures takes values close to 40%. In ad- 
dition we have found that the fraction of monomers is 
rather similar to that found for the [111] case. However, 
this is not the case for the fraction of dimers, which takes, 
in some specific cases, values below 20%. From Fig. [8] we 
also infer that there is less data dispersion for the [111] 
case than for the other two stretching directions. Note 
that part of this dispersion may be caused from the use 
of a poor statistical description for the [100] and [110] 
cases. 

When we analyze the minimum cross-section range 
0.25 < S m < 1.75 it is worth see that minimum cross- 
section histograms H(S m ) do not present a simple struc- 
ture with a single well-resolved peak centered at S m = 1 . 
On the contrary, in many cases (see Figs. we iden- 

tify the presence of satellite peaks around S m = 0.5 and 
S m = 1.5. We have explored such peaks to determine 
the type of atomic configurations included on them. In 
Fig. [5] we show the probability of finding different struc- 
tures (monomer, dimer, and complex structures) around 
each one of the sub-peaks we have commented above. For 
the sake of simplicity we have considered, for T=4K and 
T=300K, the largest nanowire sizes defined for the three 
stretching directions we have considered. In general, we 
notice that the peak at S m = 0.5 is basically originated 
from dimers whereas the peak at S m — 1 comes from 
a mix of monomer and dimer configurations. Finally, 
the region around S m = 1.5 is formed with contribu- 
tions from complex structures and a small proportion of 
monomers. This general trend is not observed at T = 4K 
for the [100] stretching direction. In this case, the differ- 



FIG. 9: Fraction of monomers, dimers and more complex 
structures that are observed around the sub-peaks located at 
S m =0.5, 1.0 and 1.5 in the H{S m ) histograms. Results for 
T = AK (left) and T = 300^ (right) are shown for the largest 
systems simulated for each one of the stretching directions 
(see Figures |3ja,b), HJa,b) and[5ja,b)). 



ence could arise from statistical errors coming from the 
short time that the system spends in the S m ~ 1 region. 

In the analysis of Figs. [8]and[9]we have argued that the 
small number of detected configurations with the Her- 
rmann " burning" algorithm is a key ingredient to under- 
stand the dispersion of values found when analyzing the 
probabilities of finding monomer, dimer or more complex 
configurations. This number depends on both the detec- 
tion window (1000 MD steps) and the intrinsic breaking 
times which could, in principle, depend on the stretch- 
ing direction, temperatures and the nanowire initial size. 
In order to get insight about the dynamics of the last 
stages of the nanowire we analyze the distribution of 
times associated to the breaking process. To do that 
we have analyzed the distribution function Po. 25 -1.75 (At) 
of the time At spent by the nanowire in the interval 
0.25 < S m < 1.75. This distribution is similar to the 
plateau lengths distribution used to experimentally de- 
tect linear atomic chains^. 

In Fig. [10] we show Po. 25-1. 75 (At) for the largest 
nanowires we used for each one of the stretching direc- 
tions as well as for three mid-size nanowires. Note that 
for the [100] and [110] cases a huge peak is observed be- 
low 5 ps. The distribution falls down fast as the time 
interval At is increased. That is, the system spends a 
short time in this minimum cross section region. For the 
[111] case we found a different dynamical pattern: the 
distribution broadens showing that the nanowire forms 
stable monomer and dimer-like configurations and giv- 
ing rise to the strong peaked H(S m ) structure depicted 
in Fig. [3] When the initial system size decreases the 
behavior for the three stretching directions is rather sim- 
ilar. In particular, the presence of two peaked structures 
in Po. 25-1. 75 (At) could correspond directly to the forma- 
tion of monomer (one atom chain) and dimers (two atoms 
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FIG. 10: Distribution function Pg, 25-1. 75 (At) of the time 
spent by the set of largest (left) and mid-size (right) Ni 
nanowires (at T = 4K) within the interval 0.25 < S m < 1.75, 
for the stretching directions [111] (bottom), [100] (middle) 
and [110] (top). 



chain), as it has been found for some metallic species^ 3 -. 
Therefore, it is important to remark that the system dy- 
namics during the last breaking stages strongly depends 
on the nanowire stretching direction and the system size. 

We have also addressed the question of the origin of 
the atoms forming the monomer and dimer configura- 
tions, in a similar way to previous studies done for gold 
nanowires^. We have performed the statistical analysis 
for all the nanowires we have considered in order to de- 
termine whether those atoms involved in the formation 
of the narrowest nanocontact region, at the last break- 
ing stages, come from bulk or surface positions. For every 
nanowire we have followed back the atomic trajectories of 
those atoms found in monomer and dimer configurations 
in order to find its initial position after the stretching 
process. We have found that surface (bulk) atoms con- 
tribute to the formation of monomers with a statistical 
weight very similar to the initial fraction of surface (bulk) 
atoms. Therefore there is not any relevant role for the 
surface atoms in the formation of monomers and dimers, 
at variance of results found for golc-2^. The origin of this 
discrepancy may arise from the comparison between two 
different metallic species (remember that, for instance, 
gold forms long atomics chains whereas nickel does not) 
or from the type and quality if the statistical analysis 
(different number of studied breaking events). 



D. Histogram decomposition into monomer and 
dimer components 

The decomposition of conductance experimental his- 
tograms into partial histograms has been already 
proposed^ 5 - to study the influence of the nanowire crys- 



talline orientation on the conductance traces or to deter- 
mine the presence of magnetic and non-magnetic contri- 
butions on Ni and Cu conductance histograms^. From a 
theoretical perspective, the decomposition on partial his- 
tograms is an additional tool that allows to understand 
the dynamical evolution of the set of nanowires taken 
into consideration to build up the full histogra m 57 ! 59 ' 73 . 

In Figs. El SI and M we plotted the minimum cross- 
section histograms H(S m ) for both temperatures and 
several nanowire sizes. We have decomposed this his- 
togram for two nanowire sizes. In black we show the 
histograms using the complete set of breaking events. In 
gray we show the partial H(S m ) histograms constructed 
from those configurations that present monomer and 
dimer configurations in the interval 0.25 < S m < 1.75. 

For the [111] stretching direction cases (see Fig. 
G2a,b,g,h)), we notice that the decomposed histograms 
are quite similar -in all the S m range- to the full his- 
togram for both temperatures and both system sizes. 
This indicates that, for the [111] stretching direction, all 
the traces S m (t) evolve from large S m values until break- 
ing using monomer and/or dimer breaking patterns. This 
breaking scheme is very robust and seems to be indepen- 
dent on the initial nanowire size and temperature (for 
temperatures below RT). 

However, for the other two stretching direction cases 
([100] and [110]) there are strong differences between 
large and small nanowire sizes. For the largest studied 
systems (see Figs. Hfa,b) and [5ja,b)), we notice that 
the monomer and dimer structures are solely responsible 
of the formation of those smaller peaks located at the 
region 0.25 < S m < 1.75 but there is a lot of config- 
urations with large S m values that do not break form- 
ing monomers neither dimers. In this case (large system 
sizes), we also see that many S m (t) curves forming long 
pentagonal nanowires (giving rise to the huge peak ob- 
served around S m ~ 5) do not break using the monomer 
or dimer breaking pattern and directly break from rel- 
atively large cross sections. The temperature does not 
play any relevant role in the studied temperatures range 
(except the rounding and smoothing of the histogram 
peaks). 

For smaller nanowires that break under [100] and [110] 
stretching (see Figs. BJg,h) and[5Ig,h)) we recover the 
breaking pattern we obtained for the [111] small size 
case, where almost all S m (t) traces present well defined 
monomer or dimer configurations during the last stages 
of the elongation process. 

Therefore, the relative formation of monomers and 
dimers during breaking events strongly depends on the 
system size when we elongate the nanowire in the [100] 
and [110] directions. The [111] stretching direction do not 
show important changes in the formation of monomer- 
dimer configurations for decreasing nanowire sizes. If 
we assume that the stretching and breaking of large size 
simulated nanowires are closer to the experimental situ- 
ation, we clearly see that the appearance of monomers 
and dimers (and the shape of H(S m )) will strongly de- 
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FIG. 11: Probability of occurrence of different atomic envi- 
ronments around a Ni monomer. A n/m configuration means 
that the monomer configuration is of type "n-l-m" or "m-1- 
n" , i.e. there are n atoms at one side of the monomer and m 
atoms at the other side. The probability has been calculated 
for the largest nanowire used for the stretching directions [111] 
(bottom), [100] (middle) and [110] (top). Two temperatures 
have been studied: 4K (black dots) and 300K (squares). 



FIG. 12: Probability of occurrence of different atomic envi- 
ronments around a Ni dimer. A n/m configuration means 
that the dimer configuration is of type "n-l-l-m" or "m-1-1- 
n", i.e. there are n atoms at one side of the dimer and m 
atoms at the other one. The probability has been calculated 
for the largest nanowire used for the stretching directions [111] 
(bottom), [100] (middle) and [110] (top). Two temperatures 
have been studied: 4K (black dots) and 300K (squares). 



pend on the experimental stretching direction used to 
elongate and break the nanowire or nanocontact. 



E. Local environment of monomer and dimer 
breaking configurations 

It is obvious that the local environment is a key factor 
to determine the conductance of monomer or dimer con- 
figurations. This has been demonstrated for Al nanocon- 
tacts where monomer and dimer transport properties 
present different feature a 57 ' 73 . For Al dimcr-like con- 
figurations, the interaction between supporting atomic 
structures surrounding the dimer becomes negligible and 
the conductance is governed by three electron channels 
whose transmission coefficients determine a conductance 
value close to Gq . 20 ^ 73 On the contrary, for a monomer 
configuration the distance between atomic supports de- 
creases, enhancing their orbital interaction, and finally 
leading to an incrase of the conductance upto G ~ 2Gn. 
We can guess a similar situation for Ni, where d orbitals 
associated to supporting pyramid-like structures will in- 
teract strongly for the monomer situation giving rise to 
higher conductance values whereas the interaction will 
decrease for the dimer-like configuration, leading to lower 
conductance values. Static conductance calculations on 
Ni configurations (monomer and dimer are consistent 
with that picture. 

From the previous analysis it is obvious that the lo- 
cal atomic environment will play an important role in 
the conductance determination, specially in the monomer 
case. MD calculations allow to establish the type of 
neighborhood formed around monomer and dimer con- 



figurations. In Figs. QT] and [12] we present the oc- 
currence probability of different atomic configurations 
around monomers and dimers. Histograms for different 
temperatures and stretching directions have very simi- 
lar shape. We can conclude that the atomic structures 
(and its probability of occurrence) around monomers 
and dimers are roughly the same independently on the 
stretching direction and the temperature. In Fig. 1111 we 
see that the most probable monomer structure presents 
the configuration 2-1-3. Around 40% of the monomers 
exhibit this configuration. But this is not the only impor- 
tant configuration since we found that 2-1-4, 3-1-3, 3-1-4 
and 2-1-2 configurations for all the stretching directions 
appear with a probability between ~ 20% and ~ 10%. 
Fig. Q2] presents the analysis of the local environment for 
those dimers we have detected in our MD simulations. 
It can be observed that the most common dimer config- 
urations are 3-1-1-4, 3-1-1-3, and 4-1-1-4 with ocurrence 
probability between ~ 40% and ~ 20%. In addition, 
some other configurations can be observed with proba- 
bility below 10%. 

It is clear that in order to analyze the experimen- 
tal conductance histograms it is important to take ac- 
count the conduction coming for all these atomic config- 
urations, including the adequate proportion of monomer 
and dimers as well as the correct proportion of possi- 
ble atomic environments for each monomer and dimer 
situation. The present information is very practical to 
address, for instance, the ab-initio study of the elec- 
tronic transport through static monomer or dimcr-like 
Ni structures. In general, the non-trivial 2-1-3 configu- 
ration would have little opportunities to be included as 
the most likely monomer structure, whereas past calcu- 
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FIG. 13: Average number of atoms < N(z) > found at a 
distance z from the monomer position, at T = 4K, for the 
stretching directions [111] (bottom), [100] (middle), and [110] 
(top). 



lations only studied 3-1-3— and 4-1- 1 29 ' 56 configurations. 
In particular, we have not found evidences of the struc- 
ture 4-1-4 in our simulations. 



F. Non-local environment of monomer and dimer 
breaking configurations 

The accurate knowledge of the local environment of 
the monomer and dimer configurations is not enough 
to understand the structure of conductance histograms. 
Electron transport mainly depends on the narrowest 
nanocontact cross-section and the surrounding regions. 
However, the electron transport is not only governed by 
the small number of open channels allowed by the nano- 
neck geometry. In fact, the transmission coefficients ap- 
pearing in Eq. ([1]) depend on several factors including 
the possible presence of disorder (defects, vacancies, im- 
purities, dislocations, ...) along the wire. 

A first approach to unveil the existence of disorder 
along the nanowire is to determine the atomic structure 
as the z distance from the monomer (or dimer) position 
increases. In Fig. ll3l we plot the average number of atoms 
< N(z) > located at a distance z from the monomer po- 
sition for the largest systems we have simulated. We have 
calculated the average number of atoms for the three 
stretching directions, at T = AK, and taking into ac- 
count all the monomer configurations appearing during 
300 breaking events. On one hand, for the [111] case, we 
note clear peaks separated a distance close to the typical 
separation of the [111] planes dm = 2.03A This points 
that the system keeps the same crystallographic struc- 
ture from the fixed slabs, at both sides of the sample, 
towards the monomer position. On the other hand, for 
the [100] and [110] stretching directions, the situation is 



quite different. Both orientations present a very similar 
local atomic structure for monomer and dimcrs. This lo- 
cal structure also agrees with that observed for the [111] 
case (as described in the previous subsection). However 
when z increases we could not distinguish a peaked struc- 
ture similat to that of the [111] case. For these cases, the 
fixed slabs try to maintain the [100] or [110] structure 
but the nanowire presents dislocations, disorder, and the 
formation of large non-crystalline wires when we study 
regions closer to the monomer. A very similar behavior 
is observed around dimmers (not shown). 

This kind of situation is not exclusive of monomer and 
dimer configurations appearing in [100] and [110] cases. 
In fact, this lack of ordering happens during the whole 
stretching process. An appealing example appears in Fig. 
[7] The unusual pentagonal structure found in the [111] 
stretching direction case (Fig. 0Jc)) presents a smooth 
transition towards the [111] supporting frozen bilayers, 
with little amount of disorder. On the contrary (Fig. 
[Zla,b)), the pentagonal configurations emerging in the 
[100] and [110] cases are connected to the frozen bilayers 
throughout rather disordered structures. 

The previously described results are very important 
whether we try to address the study of the electron trans- 
port at the latest stages of the breaking process, i.e. those 
moments when monomers and dimers govern the number 
of open channels available for electron transmission. It 
is clear that the kind of local environment of a monomer 
(or dimer) does not strongly depend on the stretching di- 
rection. Therefore the electron transport should present 
similar features for all the directions (with different oc- 
currence weights). However, the electron path inside 
the wire (neglecting inelastic collisions since the typical 
nanowire size for metals is smaller than the room tem- 
perature inelastic mean free path) could be very differ- 
ent when we compare the [111] stretching direction with 
the [100] and [110] ones. In the first case ([111]) the 
electron moves inside a well ordered system, where elas- 
tic scattering events are mainly due to the roughness of 
the nanowire walls. In the second case ([100] and [110] 
stretching directions) the electron 'feels' the presence of 
disorder during its motion in a large nanowire region. 
This disorder will likely originate a number of scattering 
events leading to a conductance decrease. If the scatter- 
ing due to disorder becomes very important the electron 
transport could leave the ballistic regime and quantum 
diffusive features should appear. Therefore we could ob- 
serve two opposite transport mechanism depending on 
the stretching direction. 

These effects possibly explain the computational con- 
ductance histograms obtained by Pauly et al^ 9 - where 
a thin nanowire (8 atoms section) was stretched along 
the [100] direction. Their calculated cross-section and 
conductance histogram shows that the contribution from 
monomer and dimer configurations in the region < 
S m < 2 is very small in comparison with the contribu- 
tions coming from the configurations S m > 2. This is 
clearly due to the chosen stretching direction (see Fig. 
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[4j. In addition, their calculated conductances also take 
values well below 3.5Go for low and relatively large cross- 
section values (with many open conducting channels due 
to the presence of s and d orbitals). This decrease of 
the expected conductance is not only due to the pres- 
ence of few atoms on the narrowest nanowire part, nor 
the d-blocking effects, but it points towards diffusive 
electron transport effects. This behavior is consistent 
with the presence of a disordered configurations along 
the nanowire, as expected for [100] stretching (see Fig. 
ITU)) . Therefore we claim that very different conductance 
histograms will obtained if thick Ni nanowires stretched 
along the [111] crystallographic direction were used. 



IV. CONCLUSIONS 

We have carried out hundreds of MD simulations of Ni 
nanowire breaking processes. We use the EAM approach 
to describe the interatomic many body interaction. From 
these simulations we are able to follow the evolution on 
time of the minimum cross-section S m (t). By adding 
hundreds of S m (t) traces we have constructed compu- 
tational minimum cross-section histograms H(S rn ) that 
statistically unveil the presence of preferred configuration 
during the elongation and breaking history. The modifi- 
cation of some initial conditions (temperature, stretching 
direction or nanowire size) leads to changes on the peaked 
structure of H(S m ) histograms. In this way, we are able 
to correlate these histogram changes with changes in the 
type of favorable structures appearing under stretching. 
The last stages of the nanowire breaking process are of 
special interest since electron transport is determined by 
a cross section formed by few atoms. However, this sit- 
uation requires an additional study since the local envi- 
ronment of those atoms forming the nanoneck plays an 
important role in the electron transport characteristics. 
Therefore, we have extracted additional information from 
the S m (t) traces as the type of configurations that appear 
at S m ~ 1. We have found that monomers, dimers and 
other more complex structures are present at the latest 
stages of the breaking events. We did not find linear 
atomic chains of three or more atoms for all the systems 
and stretching directions we analyzed. 

We have found that the H(S m ) histograms do not de- 
pend dramatically on the temperature within the ana- 
lyzed temperature range (below 300K). In general, we 
have only noticed rounding of the peaked structure of 
H(S m ), the suppression of small S m fluctuations and the 
increase, at RT, of the peak located at S m = 5 for the 
[100] and [110] cases. The absence of large thermal ef- 
fects is consistent with the fact that bulk Ni melting tem- 
perature is rather large (« 1728K), and then RT is not 
enough to activate additional mechanisms as surface dif- 
fusion effects that could modify the nanowire breaking 
dynamics. 

Summarizing our computational results we have deter- 
mined three different scenarios for the set of nanowires 



we have studied. The first scenario corresponds to large 
nanowires stretched along the [111] direction. In this 
case we have noticed a well defined series of decreasing 
H(S m ) peaks as the quantity S m increases. Nanowires 
break forming monomers and dimers, responsible of the 
large peaked structure found in the S m < 2 region. Wc 
have found that [111] breaking processes present, in gen- 
eral, well ordered structures during the breakage process, 
and monomer and dimer configurations are surrounded 
by well defined planes. This non-local ordering with re- 
spect the monomer and dimer arrangements preserves 
the ballistic character of the electron transport. This 
scenario is preserved when temperature increases from 
4K to RT. 

A second scenario corresponds to large nanowires 
stretched along [100] and [110] directions. In these cases 
the H(S m ) histogram does not large peaks in the low 
S m region (S m < 2). The analysis of such S m region 
reveals that monomer configurations also play a relevant 
in this region. However, the probability of finding dimer- 
like structures is smaller than that of finding more com- 
plex structures. In addition, the long-range environment 
surrounding monomer or dimer structures is rather disor- 
dered, possibly contributing to the conductance decrease. 
When temperature increases from 4K to RT, we found a 
dramatic increase of those peaks located at S m ~ 5 for 
both, [100] and [110] directions. We have confirmed that 
this peak is caused by the presence of long staggered pen- 
tagonal chains with ...-1-5-1-5-... structure. In relation 
with the formation of long chains in MD simulations, a 
recent work 86 points towards the higher ductibility pre- 
dicted by EAM potentials in comparison with other in- 
teratomic potentials, leading to the formation of long 
structures. However, more detailed statistical analysis 
is required to confirm this behavior. 

Finally, we can define a third scenario, common for the 
three stretching directions, corresponding to small initial 
nanowire sizes. In this case, H(S m ) histograms become 
more complicated due to size effects. This noisy scenario 
is not strongly modified when we increase temperature 
from 4K to 300K. 

In addition, we have found that monomer and dimer 
configurations appearing for the three stretching direc- 
tions present very similar features. In particular we have 
found that monomers prefer the 2-1-3 and 2-1-4 con- 
figurations to the 3-1-3 one. These configurations ap- 
pear with approximately the same probability for all the 
stretching directions and temperatures we have consid- 
ered. Most probable dimer configurations correspond to 
3-1-1-4, followed by 3-1-1-3 and 4-1-1-4 arrangements. 
These results are of importance since they define the 
kind of static configurations required to describe realistic 
breaking configurations. 

Since we have not determined the actual conductance 
associated to each atomic configuration it is very dif- 
ficult to establish a comparison between our computa- 
tional minimum cross-section histograms and experimen- 
tal conductance histograms. Furthermore, the interpre- 
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tation of experimental Ni conductance histograms is a 
difficult task since they involve not only mechanical and 
electrical properties but magnetic ones too. In spite of 
our model limitations to perform the comparison with ex- 
periments we can extract few key points from the present 
MD simulations. The similar behavior of H(S m ) at 4K 
and 300K excludes the presence of different atomic con- 
figurations as explanation of such differences between low 
and room temperatures. In addition the presence of well 
marked H(S rn ) peaks is consistent with the appearance 
of preferred conductance values (giving rise to conduc- 
tance histogram peaks). We can assume, at low tem- 
peratures and clean environment, that experimental his- 
tograms present contributions arising from breaking pro- 
cesses corresponding to different crystalline orientations. 
Thus, it is reasonable to propose that the [111] stretching 
case provides a strong contribution from monomer and 
dimer configurations, whereas [100] and [110] do not. In 
the first case (monomers and dimers) we know that ex- 
pected conductances^ fall within the range Gq~ 2.5Gq, 
just in the region where 4K experimental conductance 
histograms show a broad pea k 43 ' 46 i 51 . The second broad 
conductance peak centered around G rj 3Go could be 
explained in terms of those contributions associated con- 
figurations with S m < 2 for the three stretching diretions. 

Experimental results found at room temperature can- 
not be interpreted from the present MD simulation. The 
complex behavior of Ni conductance histograms noticed 
at room temperatures must be attributed to different 



causes: (i) the presence of structural disorder at the 
nanocontact region^, (ii) the presence of impurities or 
chemiadsorbed atoms on the nanowir o 46 ' 88 , (iii) ballistic 
magnetoresistcncc (BMR) effects originated by abrupt 
magnetic domain walls anchored at the narrowest cross- 
section of the nanowir o 35 ' 36 ! 41 ' 42 or (iv) the increase of 
the domain wall fluctuations with temperature at the 
nanocontact region. 8 ^ Therefore more realistic calcula- 
tions are required to understand the non convergent ex- 
perimental results. 
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